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Abstract 

Stochastic averaging for a class of stochastic differential equations 
(SDEs) with fractional Brownian motion, of the Hurst parameter H 
in the interval 1), is investigated. An averaged SDE for the origi- 
nal SDE is proposed, and their solutions are quantitatively compared. 
It is shown that the solution of the averaged SDE converges to that 
of the original SDE in the sense of mean square and also in proba- 
bility. It is further demonstrated that a similar averaging principle 
holds for SDEs under stochastic integral of pathwise backward and 
forward types. Two examples are presented and numerical simulations 
are carried out to illustrate the averaging principle. 
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1 Introduction 

Stochastic averaging is often used to approximate dynamical systems under 
random fluctuations. This analytic technique has been developed in the case 
of the Gaussian random fluctuations, for example, by Stratonovich [H[2] and 
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then by Khasminskii [5JH]- It has been found to be effective for understand- 
ing stochastic differential equations arising in many fields [5j [U [7] . Zhu and 
his co-workers further studied this stochastic averaging method for nonlin- 
ear systems under Poisson noise [8j [10] , and two of the present authors 
derived an averaging principle for stochastic differential equations with Levy 
noise [11] . In all of these mentioned works, the fluctuations or noises are 
uncorrelated, i.e., white noises. 

However, random fluctuations with long-range dependence, or correlated 
noises, are abundant. They may be modeled by fractional Brownian mo- 
tion (fBm) with 2 < H < 1 (where H is the Hurst index). The fractional 
Brownian motion was introduced by Kolmogorov [12]. Then, in 1968, Man- 
delbrot and Van Ness [13] presented the structure of the fractional Brownian 
motion. Due to the importance of long-range dependence of the fBm, the 
stochastic differential equations with fBm have been used as the model of the 
practical problems in various fields, such as hydrology, queueing theory and 
mathematical finance (Chakravarti and Sebastian, p3]; Hu and 0ksendal, 
[To] ; Leland, Taqqu, Willinger, and Wilson et al, [16]; Scheffer, [17]). So 
fractional Brownian motion has also been suggested as a replacement of 
standard Brownian motion in several stochastic models (|18 [ \19 \ 120]). 

Given the abundance of correlated fluctuations, it is crucial to un- 
derstand the behaviors of the stochastic differential equations with fBm 
|2H [22| [231 [23] . Unfortunately, the fractional Brownian motion is neither a 
semi-martingale nor a Markov process, so the powerful tools for the stochas- 
tic integral theories are not applicable when studying fBm. Therefore, much 
of the recent research on SDEs with fBm is by numerical simulations. Other 
techniques for such SDEs would be desirable. This motivates us to inves- 
tigate stochastic averaging techniques for differential equations driven by 
fractional Brownian motion. 

In the present paper, we study a stochastic averaging technique for a 
class of SDEs with fBm. We present an averaging principle, and prove 
that the original stochastic differential equation can be approximated by 
an averaged stochastic differential equation in the sense of mean square 
convergence and convergence in probability, when a scaling parameter tends 
to zero. In addition, the similar conclusion holds for a SDE, where the 
stochastic differential or stochastic integral is of forward and backward types. 

The organization of the paper is as follows. Section 2 recalls the defi- 
nition of fractional Brownian motion and highlight the differences with the 
usual Brownian motion roughly, and then briefly reviews the symmetric, 
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forward and backward stochastic integrals with respect to fBm. Section 3 
is devoted to prove a stochastic averaging principle for stochastic differen- 
tial equations with fBm. Section 4 presents two examples to illustrate the 
stochastic averaging principle. 

2 Fractional Brownian motion and stochastic in- 
tegration 

Since stochastic differential equations are interpreted via stochastic inte- 
grals, it is necessary to specify the integration with respect to fBm. For 
background on this issue, see [251 [23 123 123 122] • For instance, using the 
notions of fractional integral and derivative, it is appropriate to introduce a 
pathwise stochastic integral with respect to fBm [6U \ \'61\ 152]. 

In this preliminary section, we briefly recall the definition of fBm and 
the integration with respect to it, for H £ (|, 1). 

2.1 Fractional Brownian motion 

Let (f2, J£", P) be a complete probability space. The definition of the frac- 
tional Brownian motion is as follows [13"]. 

Definition 1. The fractional Brownian motion (B H (t)) with Hurst index 
H is a centered self-similar Gaussian process B H = B H (t),t£ R+, on 
(0, P) with the properties : 

(1) B H (0) = ; 

(2) EB H (t) =0,t£K + ; 

(3) KB H {t)B H {s) = ±(\t\ 2H + \s\ 2H -\t- s\ 2H ),t, s€R + . 
For H = h, this is the usual Brownian motion. 

We also recall the following features of the fractional Brownian motion: 

(a) Self-similarity : For every constant a > and every T > 0, the 
following relation about distribution (or law) holds 

Law {B H {at) : t G [0, T]) = Law (a H B H (t) :te[0, T]). 

The above formula means that the two processes B H (at) and a H B H (t) have 
the same finite-dimensional distribution functions, i.e., for every choice of 
to, . . . ,t n £ 

P(B H (at ) <x ,..., B H {at n ) < x n ) = P(a H B H (t ) <x ,..., a H B H (t n ) < x n ), 
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for every xo, . . . , x n € R. 

(b) Stationary increments : The increment of this process in (s, t) has a 
normal distribution with zero mean, and the following variance 

E{B H (t) - B H (s)) 2 = \t-s\ 2H . 

Hence, for every integer k > 1, 

E(B H (t)-B H (s)f k = ^\t-s\ 2H K 

In other words, the parameter H controls the regularity of the trajectories. 
For H = 2, the increments of the process in disjoint intervals are indepen- 
dent, while for H ^ ^, the increments are dependent. 

(c) Long-range dependence : The auto-covariance function pH(n),n £ N 
of the fBm is 



p H (n) ■= Cov(B H (k) - B H (k - 1), B H (k + n) - B H (k + n - 1)) 

2n 2H ] 

H(2H - l)n 2H - 2 , 



\(n+l) 2H + (n-l) 2H -2n 2H ] 



and p_ff(n) — > 0, as n tends to infinity. 

If H > 2, Ph(ti) > 0, for n large enough, and Y^=i PH{n) = oo. In 
this case, we say that the fractional Brownian motion has long-range depen- 
dence. So the fBm can be used to describe cluster phenomena, occuring in 
geophysics, hydrology and economics. 

Based on the definition of the fractional Brownian motion, it is clear that 
the standard Brownian motion is a specific fractional Brownian motion with 
index H =1/2. 

The relationship between the usual Brownian motion and fractional Brow- 
nian motion is as follows: 

(Rl) The similarities : They are both Gaussian process; they do not have 
differentiable sample paths and both have statistical self-similarity; besides 
they are almost everywhere Holder continuous. 

(R2) The differences : Fractional Brownian motion is neither a semi-martingale 
nor a Markov process (for H ^ i), but the usual Brownian motion is a 
semi-martingale and a Markov process; fractional Brownian motion has no 
independent increments, while the usual Brownian motion has. 
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2.2 Stochastic integration with respect to fractional Brown- 
ian motion 

For the convenience of readers, we recall some stochastic integration with 
respect to the fractional Brownian motion [33j G21 E5] • 

Let <p : R + x R + — > E + be given by 

ip(t, s) = H(2H — — s\ 2H ~ 2 , t, s £ R+, 

where i < H < 1, and let / : R + — > M + be Borel measurable. Define 
L% = {f:\f\l= I [ f(t)f(sMt,s)dsdt<oo}. 

The Hilbert space 1?^ is naturally associated with the Gaussian process 
(B H (t),t>0). 

Let S be the set of smooth and cylindrical random variables of the form 

F = f(B H (ip 1 ),B H & 2 ),...,B H W n )) ) 

where n > 1, / € C£°(lR n ) (i.e., / and all its partial derivatives are bounded), 
and ipi £ H, H is a Hilbert space [29] . 
Introduce the Malliavin ip — derivative of F 

DfF = [ <p(t, u)D H Fdv, 

where 

i=l 1 

In this paper, we consider the pathwise stochastic integrals for fBm. The 
definition of the symmetric stochastic integral for the fBm case is in [33]. 

Definition 2. Let u(t) (t £ [0, T]) be a stochastic process with integrable 
trajectories. The symmetric integral of u(t) with respect to B H (t) is defined 
as 

limj- f T u(s)[B H (s + e)-B H (s-e)}ds, 
provided that the limit exists in probability, and is denoted by J Q T u(s)d° B H (s) . 
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Remark 1. Let L(0,T) be the family of processes on [0,T], such that u(t) £ 
L(0,T) j/E|tt(i)|?, < oo. Assume that u(t) (t € [0, T]) is a stochastic process 
in L(0, T) and satisfies 



[ [ \Df u(t)\\t - s\ 2H - 2 dsdt < oo. 
Jo Jo 

Then the symmetric integral exists and the following relation holds: 

u(t)d°B H (t)= [ u(t)odB H (t) + [ Dfu(s)ds, 
Jo Jo 



where o denotes the Wick product, H G (^,1)- 

Remark 2. The definition of the forward and backward integrals with respect 
to fBm is as follows: 

Let u(t) (t G [0, T]) be a process with integrable trajectories. The forward 
integral of u(t) with respect to B H (t) is defined as 

e^o e Jo e 

provided that the limit exists in probability, and is denoted by J Q T u(s)d~ B H \s) . 
The backward integral is defined as 

1 f T ^ B H (s-e)-B H (s) ^ 

lim - / u(s)\ \ds, 

e^o e J e 

provided that the limit exists in probability, and is denoted by J T u(s)d + B H (s) . 



Remark 3. According to \33f . under the assumptions in Remark 1, the 
symmetric, backward and forward integrals coincide in the following sense 



[ u(t)d~B H (t)= [ u(t)odB H (t)+ f Dfu(s)ds, (2) 

Jo Jo Jo 

[ u(t)d + B H (t)= [ u(t)odB H (t)+ f Dfu(s)ds. (3) 

Jo Jo Jo 



I. 



3 An averaging principle for SDEs with fBm 

3.1 Some Lemmas 

In order to present a stochastic averaging principle, we need two lemmas. 



Lemma 1. Let B H (t) be a fractional Brownian motion with \ < H < 1, 
and Z(s) be a stochastic process in L[0,T] . For every T < oo , there exists 
a constant C(H,T) = jjx 2 ^ 1-1 such that 

E[( f T \Z{s)\ o dB H {s)f] < C(H, T)E[ f \Z(s)\ 2 ds] + CT 2 . (4) 
Jo Jo 

Proof. According to [33] ( Theorem 2.1 ), 



E f (Df\Z(s)\) 2 ds < oo, 
Jo 



and 



E( f T \Z(s)\ o dB H (s)) 2 = E[\\Z\\l + f Df\Z( S )\ds) 2 ]. 
Jo Jo 

Thus, 

E([ \Z(s)\odB H (s)) 2 = E[ I [ Z(t)Z(s)(p(t,s)dsdt + M(f DfZ(s)dsf 
Jo Jo Jo Jo 



^A + B, 

where 

A = e£ fT\Z(t)\\Z( S )\ip(t,s)dsdt, 
B = E(tf Df\Z(s)\ds) 2 . 

We further have 

T r T 

2H ~ 2 dsdt 



A < H(2H - 1)E f [ \Z(s)\ 2 \t 
Jo Jo 

<H(2H-1)E^ |Z( S )| 2 ^-_ T 



ds 



T 



< HT 211 '^ / \Z(s)\ 2 ds 
Jo 

r T 

<C(H,T)E \Z{s)\ 2 ds. 
Jo 
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By the Cauchy-Schwarz inequality for B , we get : 

rT 



B<T 2 e[ \Df\Z(s)\\ 2 ds <CT 2 
Jo 
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Then we can finally deduce that 

E[( [ T \Z(s)\odB H (s)) 2 } < C(H,T)E C \Z \s)\ 2 'ds + CT 2 . 
Jo Jo 

This finishes the proof of this Lemma. 

Lemma 2 can be obtained according to Definition 4 and Lemma 1. 

Lemma 2. Suppose that Z(s) is a stochastic process in L[0, T], and B H (t)(H > 
|) is a fractional Brownian motion. For any < T < oo, there exists a con- 
stant C{H,T), such that the following inequality holds 

E[( [ T \Z(s)\d°B H (s)) 2 } < 2C(H, T)E[ C \Z(s)\ 2 ds] + ACT 2 , (5) 
Jo Jo 

where C(H, T) = HT 2H ^ . 

Proof. Using Eq.(l) and the Cauchy-Schwarz inequality, we can get 

E[( f \Z(s)\d°B H (s)) 2 ] = E[( [ T \Z(s)\ o dB H {s) + f Df\Z(s)\ds) 2 ] 
Jo Jo Jo 

< E[2( f \Z{s)\ o ^(,)) 2 + 2( [ T Df\Z(s)\ds) 2 ] 
Jo Jo 

<2E[(/ \Z{s)\ o B H (s)) 2 } + 2E[ f Df\Z(s)\ds% 
Jo Jo 

Due to Eq.(4) and 

E( / Df \Z{s)\ds) 2 < CT 2 , 
Jo 

we obtain 

E[(/ \Z(s)\d°B H (s)) 2 } < 2C(H,T)E / \Z{s)\ 2 ds + 2CT 2 + 2CT 2 , 
Jo Jo 

namely 

E[( [ T \Z(s)\d°B H (s)) 2 } < 2C(H, T)E[ C \Z(s)\ 2 ds] + ACT 2 . 
Jo Jo 

the proof is completed. 



3.2 Stochastic differential equations driven by fractional Brow- 
nian motion 

In this section, we concern the symmetric integral of stochastic differential 
equations with respect to fBm. Solutions of the stochastic differential equa- 
tion driven by fractional Browinan motion have been studied intensively by 
using the pathwise approach [36j EZ] • 
Consider the equation on R d 

X(t)=X(0)+ [ b{s,X(s))ds+ [ a(s,X(s))d°B H (s), (6) 
Jo Jo 

where X(0) is a given d-dimensional random variable, b(s,X(s)) : [0, T] x 
M. d — > Mr is a measurable vector function, a(s, X(s))isad x m matrix with 
each element aj t i(s, X(s)) : [0, T] x M d — > M d a measurable vector func- 
tion, and the processes B H (t), represents d-dimensional fractional Brown- 
ian motions with Hurst parameter H defined in a complete probability space 
(O, J 7 , P). Denote by a = {c?j : i) d '™ =l the matrix of "diffusion" and b = (bi)f =1 

the "drift" vector, \a\ := (£■ ^ \a hl \ 2 )^ , \b\ := (E,(&i) 2 )^ • 

Let us consider the following assumptions on the coefficients : 

a(t, x) is differentiable in x , and satisfies : there exists M > , < 7 , 
k < 1 , and for any N > 0, M N > 0, 

(i) a is Lipschitz continous in x ,V x, y G M. d , t £ [0, T] : 

\a(t,x) - a(t,y)\ < M\x - y\, 

(ii) x -derivative of a is local Holder continous in x , V |y| < N, t 6 [0, T] 

\a x (t,x) - a x (t,y)\ < M N \x - y\ k , 

(iii) <7 is Holder continous in times, for all x G M. d ,t, s G [0, T] : 

\a(t,x) - cr(s,x)\ + \a Xi (t,x) - a Xi (s,x)\ < M\t - s| 7 . 
for each i = 0, . . . , d. 

The function b = b(t, x) satisfies the following conditions: 

(iv) for all N > , there exist L N > , for all \x\,\y\ < N,t G [0,T] , such 
that 

\b(t,x) - b(t,y)\ < L N \x - y\, 
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(v) there exists the function 60 £ L p (0,T;M. d ) (p > 2) , and L > , for any 
(t, x) £ [0, T]xi such that 

\b(t,x)\ < L\x\ + b (t). 

On the basis of Theorem 3.1.4 and Remark 3.1.5 in [35], there exists the 
unique solution (X t ,t £ [0,T]) of the Eq.(6). 

3.3 An averaging principle 

Now we discuss a standard stochastic differential equation using an averag- 
ing principle in M. d . 

The standard stochastic differential equation is defined as: 

X € (t) = X(0) + e 2H f b{s,X t {s))ds + e H f a(s, X e (s))d° B H \s). (7) 
Jo Jo 

where X(0) = Xq is a given d-dimensional random varibale as the initial 

condition, t £ [0, T] and the coefficients have the same conditions as in 

Eq.(6), and e £ (0, eo] is a positive small parameter with eo a fixed number. 

Assume that (i) — (v) (the Lipschitz and growth conditions) are satisfied, 

besides the mappings b : M. d — > M. d , a : M. d — > M. d , are measurable. And 

presume they meet the following additional inequalities : 

(CI) 

^ [ Tl \b(s,y) - b(y)\ds < MTi)(l + |y|), 
J i Jo 

(C2) 

^ [ Tl \a(s,y) - a(y)\ 2 ds < ^(Ti)(l + \y\ 2 ). 
J i Jo 

where T\ £ [0, T], fi(Ti) are positive bounded functions with lim^^oo (fi(Ti) 
= 0, i = 1,2. 

Then, we can obtain the SDEs with the averaging principle : 

Z e (t) = X(0) + e 2H f b(Z € (s))ds + e H f a(Z e (s))d° B H \s). (8) 
Jo Jo 

This SDE is called the averaged SDE of the original standard SDE (7). Un- 
der the similar conditions such as X(t) in Eq.(6), this equation will have a 
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unique solution Z e (t). 

Now We claim the following main theorems to show relationship between 
solution processes X e (t) and Z e (t). It shows that the solution of averaged 

Eq.(8) converges to that of the original Eq.(7) in the sense of mean square 
and probability respectively. 

Theorem 1. Suppose that the original SDEs (7) and the averaged SDEs (8) 
both satisfy the assumptions (i)-(v)and (C1)-(C2). For a given arbitrarily 
small number 8\ > , there exist L > , e\ G (0, eo] and j3 G (0, 1), such 
that for any e G (0, ei] , 

E(\X e (t) - Z t (t)\ 2 ) < Si- 

Remark 4. (i) This conclusion shows that the solution of averaged SDEs 
converges to that of initial SDEs in a certain sense. That is Theorem 1 
means the convergence of these two solutions in the sense of mean square. 

(ii) If only partial conditions hold, Theorem 1 may still hold. In this 
situation we may speak of partial averaging. 

Proof. 

According to the above analysis, we start with 

X e {t) - Z e (t) = e 2H [ [b(s,X e (s)) - b(Z t (s))]ds + 
J o 

e H f\a(s,X e (s))-a(Z e (s))]d°B H (s), 
Jo 

and employ the following inequality for n G N, and x±,X2, ■ ■ ■ ,x n G R : 

\xi + x 2 + ■■■ + x m \ 2 < m(|xi| 2 + |x 2 | 2 + • • • + l^ml 2 ), (9) 
we arrive at 

\X&) - Z e {t)\ 2 < 2e 4H \ [\b(s, X e (s)) - b(Z € (s))]ds\ 2 + 

Jo 

2e 2H \ [\a(s,X e ( S ))-HZe(s))]d°B H ( S )\ 2 
Jo 

= /? + /*. 

where [0, t] G [Q,u] G [0, T],Ii,i = 1,2 denote the above terms respectively. 
Now we present some estimates for Ii, i = 1,2. 
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Firstly, we apply the inequality (9) to get 

I 2 = 2e m \ [ [b(s,X e (s))-b(Z £ (s))]ds 
Jo 



ft 

AH 



<2e 4 "| / [b(s,X e (s))-b(s,Z e (s)) + b(s,Z £ (s))-b(Z e (s))]ds\ 



o 



<4e 4H | / [b(s,X e {s))-b(s,Z £ {s))]ds\ 2 + 
Jo 

Ae 4H \ f[b(s,Z e {s))-b(Z e (s))]ds\ 2 
Jo 



where 



I 2 n = 4e 4H \ [ [b(s,X e (s))-b( S ,Z e (s))]ds\ 2 , 
Jo 

I 2 2 = Ae 4H \ f\b(s,Z e (s)) - b(Z e (s))]ds\ 2 . 
Jo 

By the Cauchy-Schwarz inequality for I 2 ^ , we obtain : 

|/n| 2 < ±e AH t\ f\b(s,X £ (s)) - b(s,Z e (s))} 2 ds\, 
Jo 

Because of condition (ii) and taking expectation, we can get 

E|/n| 2 < Ae 4H E(t f \b{s, X e (s)) - b(s, Z e {s))\ 2 ds) 
Jo 

<4e 4H uL 2 N E [ \X e (s) - Z e (s)\ 2 ds) 
Jo 

< 4e 4H uL 2 N E( [ U \X £ (s) - Z e (s)\ 2 ds) 

Jo 

i-u 

< Ae AH uKu / E(|X e (si) - Z e (si)\ 2 )du. 

Jo 

where K\\ is a constant. 

Then about l\ 2 , we use condition (CI), (fi(t) is positive bounded func- 
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tion and take expectation to yield : 

E|/i2| 2 < 4e m E(t 2 [- f \b(s, Z e (s)) - b(Z e (s))\ds} 2 ) 
t Jo 

<8e AH u 2 E{ Vl {t) 2 (l + \Z e {s)\ 2 ))} 
<8e 4H u 2 K 12 (l+E( sup (\Z e (t)\ 2 )) 

0<t<u 

< 8e* H u 2 K 12 , 

where K\ 2 denotes a constant which may differ in the above inequality. For 
each t > , we get 

E|/i| 2 < 4e 4H uK n / E(\X t ( Sl ) - Z e ( Sl )\ 2 )du + 8e m u 2 K l2 . {Z x ) 
Jo 

Now take expectation on l| to obtain 

E|/ 2 | 2 = 2e 2H E\ f\a(s,X e ( S )) - a(Z € (s))]d° B H (s)) 2 \ 
Jo 

< 4e 2H E(\ f\a(s, X € (s)) - a(s, Z € (s))]d° B H (s)\ 2 ) + 
Jo 

4e 2H E(\ [\a(s,Z e ( S ))-a(Z e (s))]d B H (s)\ 2 ) 
Jo 

_ t2 _|_7-2 
— J 21 + i 22- 

where 

1 2 21 =4e 2H E(\ [\a(s,X t (s)) - a(s,Z e (s))]d°B H (s)\ 2 ), 

Jo 

1 2 22 = 4e 2H E(\ [\a(s, Z e (s)) - a(Z e (s))]d° B H (s)\ 2 ). 

Jo 

By the Lemma 2, conditions (i) and (C2), it is easy to get 

II < 4e 2H E(\ f M\X e (s) - Z e (s)\d°B H ( S )\ 2 ) 
Jo 

<4e 2H M 2 E\ [ \X e (s)- Z e (s)\d°B H (s)\ 2 
Jo 

< 4e 2H M 2 {2Ht 2H ~ l E[ [ \X e (s) - Z e (s)\ 2 ds] + 4Ct 2 ) 

Jo 

< 8e 2H Hu 2H ~ l K 21l [ E{\X t ( Sl ) - Z t ( Sl )\ 2 )du + 16e 2H u 2 K 212 . 

Jo 
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Due to the conditions (C2) , we obtain 

I 2 22 < 4^E(| [ U [a(s,Z e (s))-a(Ze(s))]d°B H (s)\ 2 ) 
Jo 

< Ae 2H {2Hu 2H ^E[ [ \a(s, Z € (s) - a(Z e (s))\ 2 ds] + ACu 2 } 

Jo 

< B^Hu 211 -^ [ \a(s, Z e {s) - a(Z £ (s))\ 2 ds] + 16e 2H Cu 2 

Jo 

< 8e 2H Hu 2H K 22l {[l+E( sup \Z e (t)\ 2 )}} + 16e 2H Cu 2 

0<t<u 

< 8e 2H u 2H HK 221 + 16e 2H u 2 K 222 . 

where the last inequality is obtained by the same arguments of E|7i| 2 , and 
K 2 i,i = 1,2 denote positive constants that may differ in different cases. 
Then 

E|I 2 | 2 < Ae 2H Hu 2H ~ l K 211 [ E(\X e ( Sl ) - Z t ( Sl )\ 2 )du + 

Jo 

\Qe 2H u 2 K 2l2 + 8e 2H u 2H HK 22l + 16e 2H u 2 K 222 . (Z 2 ) 
Therefore from above discussions (Z±) and (Z 2 ), we can get 
E(\X € (t) - Z e (t)\ 2 ) 

< 4e 4H uK u [ R(\X £ (si) - Z e { Sl )\ 2 )du + 8e m u 2 K 12 + 

Jo 

8e 2H Hu 2H ~ l K 2ll / E(|X e (si) - Z e ( Sl )\ 2 )du + We 2H u 2 K 212 + 
Jo 

8e 2H u 2H HK 221 + 16e 2H u 2 K 222 

< (Ae^uKn + Se^Hu 2 *- 1 ^) f ' M(\X e ( Sl ) - Z t ( Sl )\ 2 )du + 

Jo 

8e AH u 2 K l2 + 16e 2H u 2 K 2l2 + 8e 2H u 2H HK 22l + l6e 2H u 2 K 222 . 
Now by the Gronwall-Bellman inequality, we obtain 
E(|X e (t) - Z e (t)\ 2 ) 

< 8e 2H u(e 2 uK 12 + 2uK 212 + u 2H ~ l HK 221 + 2uK 222 ) 
exp4e 2H (e 2 uK n + 2Hu 2H ~ 1 K 211 ). 



Select p G (0,1), L > , such that for all t G (0,Le~ 2H/3 ] C [0, T] , we 
have 

E(|X e (i) - Z e (i)| 2 ) < ifae 1 - 2 ^. 
14 



where 

K 3 = 8Le l - 2H P(Le 3 - 2H ?- 2H K 12 + 2Le l ~ 2H P- 2H K 212 + 

L 2H-l e 2Hp-W^-AH* + AH-l HK22i + 2L e 1 - 2H ^ 2H K 222 ) 

exp4e 2H (L e 3 - 2 ^- 2 ^ n + ^-1^-4^/3-2^+2^+4^-1^) 
it is a constant. 

Consequently, given any number b\ > , we can select e\ £ (0, eo] , such 
that for every e £ (0, ei] , and for each t E (0, Le~ 2H/3 ] 

E(\X e (t) - Z e (t)\ 2 ) <Si. 

This is all of the proof. 

We also have the following result on uniform convergence in probability. 

Theorem 2. Suppose that all assumptions (i)-(ii) and (C1)-(C2) are sat- 
isfied. Then for any number 5 2 > , we have 

IimP(|X e (t)-Z e (t)| >5 2 ) = 0, 

where L and are the same to Theorem 1. 
Proof. 

On the basis of Theorem 1 and the Chebyshev-Markov inequality, for 
any given number 5 2 > , one can find 

F(\X € (t) - Z e (t)\ >5 2 ) 

<^E(\X e (t)-Z e (t)\ 2 )<^ 2H "- 
Let e — > and the required result follows. 

Remark 5. Theorem 2 means the convergence in probability between the 
original solution X e {t) and the averaged solution Z t {t) . 

Then, we also can study the forward integral and backward integral of 
stochastic differential equations driven by fBm, and the definition of the 
forward integral and backward integral are the same to section 2 : 

dX(t) = b(t, X{t))dt + a(t, X(t))d-B H (t), (10) 

dX(t) = b(t, X{t))dt + a(t, X{t))d + B H (t). (11) 
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On the basis of the Eq.(7) and Eq.(8), we can get the standard stochastic 
differential equation and the averaged SDEs : 



X £ (t) = X{0)+e 2H f b{s,X e {s))ds+e H f a{s, X e (s))d~ B H (s), 
Jo Jo 



(12a) 



Z e (t) =X(0)+e 2H [ b(Z t (s))ds + e H [ a(Z e (s))d~ B H (s), (126) 
Jo Jo 

X t {t) = X(0)+e 2H [ b{s,X e {s))ds+e H [ a(s, X e {s))d + B H (s), (13a) 
Jo Jo 

Z e {t) = X(0) + e 2H f b{Z e {s))ds + e H [ a{Z € (s))d + B H \s). (136) 
Jo Jo 

where X(0) = Xq is the initial condition, and the coefficients satisfy the 

{%) — (v) conditions. 



Theorem 3. Assume the original SDEs (12, 14) and the averaged SDEs 
(13,15) both satisfy the (i)-(v) and (C1)-(C2). For a given arbitrarily small 
number 63 > , there exists Q > , €2 £ (0, eo] and f3 £ (0,1) , such that 
for any e G (0, e 2 ] , 

E(\X € (t) - Z e (t)\ 2 ) < 
And then for any number 84 > 0, we can get 

\im¥(\X e (t) - Z e {t)\ > 5 4 ) = 0. 

e— >0 

Proof. 

Due to the Theorem 1 and Theorem 2, this proof is the similar to the process 

of SDEs with the symmetric integral. 

We regard the forward integral of SDEs as an example. 

X e (t) - Z e (t) = e 2H [ [b(s,X e (s)) - b(Z e (s))]ds + 
Jo 

e 11 [ t [a(s,X e (s))-a(Ze(s))}d-B H (s), 
Jo 



Then we can obtain 



X e (t)-Z e {t)\ 2 <2e AH \ [ [b(s,X e (s))-b(Z e (s))]d S \ 2 + 

Jo 

2e 2H \ f[o-( Sl X e (s))-a(Z e ( S ))]d~B H { S ^ 
Jo 



F 2 + F 2 , 
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where < t < u G [0, T],Fi, i = 1, 2 denote the above terms respectively, we 
get 

Ft <4e 4H \ [ t [b(s,X e ( S ))-b(s,Z e (s))}ds\ 2 + 
Jo 



Ae 4H \ I [b(s,Z e (s))-b(Z e ( S ))]ds\ 



o 

<F 2 1+ F 2 2 , 
The similar technique yields 



l^iil 2 <4e 4H t\ [\b(s,X € ( S ))-b(s,Z € (s))] 2 ds\, 
Jo 

E|Fn| 2 < 4e 4// E(t [ \b(s,X e {s)) - b(s, Z e (s))\ 2 ds) 
Jo 

< Ae 4H uL u [" E(\X e ( Sl ) - Z € ( Sl )\ 2 )du, 
Jo 



where L\\ is a constant. 

r-t 



E\F l2 \ 2 <Ae iH E(t 2 [\ I \b(s,Z e (s))-b(Z e (s))\ds] 2 ) 
t Jo 

< 8e AH u 2 L 12 , 

where L\ 2 denotes a constant which may differ from L\\. 
We obtain 

IE I i 7 ! 1 2 < Ae 4H uL\\ f\{\X t { Sl )- Z e { Sl )\ 2 )du + %e AH u 2 L l2 . (Fx) 
Jo 

Consider the F 2 , 

E\F 2 \ 2 < Ae 2H E(\ [\a(s,X e (s))-a(s,Z e ( S ))}d°B H ( S )\ 2 ) + 
Jo 

Ae 2H E(\ f\a(s, Z e {s)) - a(Z e (s))]d° B H (s)\ 2 ) 
Jo 

= fI x + f| 2 , 

By the previous conditions, we can get 

F 2 X < Ae 2H E{\ I M\X e {s) - Z e (s)\d~ B H \s)\ 2 ) 
Jo 

< Se 2H Hu 2H - l L 2ll / E(\X t { Sl ) - Z t { Sl )\ 2 )du + 16e 2H u 2 L 212 , 
Jo 
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F| 2 < 4e 2// E(| [ U [a(s, Z e (s)) - a(Z e (s))]dT B H \s)\ 2 ) 
Jo 



< 8e 2H u 2H HL 221 + \Qe 2H u 2 L 222 . 

Therefore 

E|F 2 | 2 < Ae 2H Hu 2H ~ l L 21l / E(|X e (*i) - Z £ ( Sl )| 2 )cb + 

lQe 2H u 2 L 212 + 8e 2H u 2H HL 221 + 16e 2// n 2 L 222 . (y 2 ) 
Considering (Yi) and (Y 2 ), one arrives at 

E(\X e (t) - Z e (t)\ 2 ) < (^ H uL u + %e 2H Hu 2H - l L 2ll ) / E(|X e ( Sl ) - Z e { Sl )\ 2 )du + 

JO 

8e 4// u 2 L 12 + 16e 2H u 2 L 2l2 + 8e 2H u 2H HL 221 + 16e 2H u 2 L 222 . 

The discussions that follow are same to the process of proofs to Theorem 
1 and Theorem 2. 

Remark 6. That is to say, by Theorem 1, Theorem 2 and Theorem 3, we 
can get the same results for three types of pathwise integrals of SDEs. 

4 Examples 

Through the above discussion, we have established an averaging principle 
for the SDEs (6) with fractional Brownian motion. For Eq.(7) we can define 
the standard SDEs and the averaged SDEs respectively 

dX e (t) = e 2H b(t,X € (t))dt + e H a(t,X,(t))d°B H (t), (12) 

dZ e {t) = e 2H b(Z e (t))dt + e H a(Z e (t))d°B H (t). (13) 
with the same initial condition 

X e {0) =Z e (0) = x . 

Assume that the conditions of Theorem 1 are satisfied for b, a , and 
the similar conditions (C1)-(C2) are satisfied for b, a. Then the following 
averaging principle holds 

K(\X e (t) - Z e (t)\ 2 ) <5 U 
limP(|X £ (t) - Z € (t)\ >8 2 ) = 0. 

e->0 
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where the constants L, f3, e, 6%, 82 are the same as in Theorem 1 and The- 
orem 2. 

Now we present two examples to demonstrate the procedure of the av- 
eraging principle. 

Example 1. Consider the following SDEs driven by fractional Brownian 
motion : 

dX e = -2e 2H \X e sin 2 (t)dt + e H d°B H {t), (14) 

with initial condition X e (0) = Xo and E,\Xo\ 2 < 00, where b(t,X e ) = 
—2XX e sm 2 (t),a(t,X e ) = 1, and X is a positive constant, B H (t) is a frac- 
tional Brownian motion. Then 

b(X e ) = - b(t, X e )dt = -\\X e , a(X € ) = 1, 
^ Jo 2 

and define a new averaged SDE 

dZ e = e 2H b(Z e )dt + e H a(Z e )d°B H {t), 

namely, 

dZ € = -\? u \Z e dt + e H d°B H (t), (15) 

Obviously, Z e (t) is the well-known Ornstein-Uhlenbeck process, and the so- 
lution can be obtained as : 

Z e (t) = exp(-ieAt)X + yfe [ exp(-i e A(t - s))d°B H {t). 
2 Jo 2 

Because all the conditions (i)-(ii) and (C1)-(C2) are satisfied for func- 
tions b,a,b,a in SDEs (6), (7), thus Theorem 1 and Theorem 2 hold . 
That is, 

E(\X e (t) - Z e (t)\ 2 ) <«5i, 

and as e — > , 

X e (t) — > Z e (t) in probability. 

Now we carry out the numerical simulation to get the solutions of (|14p 
and (US]) under conditions of X = 0.0, A = 0.2, e = 0.045, H = 0.75 , 
X = 0.1, A = 0.2, e = 0.045, H = 0.85 ,X = 0.1, A = 0.4, e = 0.01, H = 0.6, 
X = 0.0, A = 0.4, e = 0.02, # = 0.7 respectively. Figure 1 (a) ~ (d) 
show the comparison of exact solution X t (t) with averaged solution Z t (t) for 
equations (|14p and (|15p . One can find a good agreement between solutions 
of the original equation and the averaged equation. 
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Figure 1: Comparison of the exact solution X e (t) with the averaged solu- 
tion Z t (t) for equations Eq.(14) and Eq.(15) (a) X = 0.0, A = 0.2, e = 
0.045, H = 0.75 (b) X = 0.1, A = 0.2, e = 0.045, H = 0.55 (c) X = 0.1, A = 
0.4, e = 0.01, H = 0.6 (d) X = 0.0, A = 0.4, e = 0.02, H = 0.7 



Example 2. Consider the following SDE with the fractional Gaussian noise 

dX e = -e 2H dt + e H cos 2 (t)\d°B H {t), (16) 

here we denote X(0) = Xq as the initial condition with E|X(0)| 2 < oo. 
Here b(t,X e ) = —1, and a(t,X e ) = cos 2 (i)A. Now we define a new (av- 
eraged) SDE as 

dZ e = -e 2H dt + ^Xe H d°B H (t). (17) 

Obviously, all conditions in Theorem 1 and Theorem 2 are satisfied 
for the averaged SDE (17), so we can use the solution Z e (t) to approximate 
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the original solution X e {t) to SDE (16), and the convergence in mean square 
and in probability will be assured. 
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Figure 2: Comparison of the exact solution X £ (t) with the averaged solution 
Z £ (t) for equations(17).(a) X = 0.0, A = 2.0, e = 0.001, if = 0.55 (b) X = 
0.0, A = 2.0,e = 0.0045, H = 0.65 (c)X = 0.1, A = 3.0, e = 0.002, if = 0.6 
(d).Xo = 0.0, A = 3.0, e = 0.002, H = 0.7 



Acknowledgments 

This work was supported by the NSF of China (Grant Nos. 10972181, 
11102157), Program for NCET, the Shaanxi Project for Young New Star in 
Science and Technology, NPU Foundation for Fundamental Research and 
SRF for ROCS, SEM. We thank Ilya Pavlyukevich for valuable discussions, 
and referees for their helpful comments. 



21 



References 



[I] R. L. Stratonovich, Topics in the Theory of Random Noise, New York, 
Gordon and Breach, Vol.1, 1963; Vol.2, 1967. 

[2] R. L. Stratonovich, Conditional Markov Processes and Their Application 
to the Theory of Optimal Control. American Elsevier, 1967. 

[3] R. Z. Khasminskii, A limit theorem for the solution of differential equa- 
tions with random right-hand sides, Theory Probab.Appl. 11(1963), 390- 
405. 

[4] R. Z. Khasminskii, Principle of averaging of parabolic and elliptic 
differential equations for Markov process with small diffusion, Theory 
Probab.Appl.8(l)(1963) 1-21. 

[5] N. Sri. Namachchivaya and Y. K. Lin , Application of stochastic averag- 
ing for systems with high damping, Probab. Eng. Mech. 3(1988) 185-196. 

[6] J. Roberts and P. Spanos, Stochastic averaging : an approximate method 
of solving random vibration problems, Int. J. Non- linear Mech. 21(1986) 
111-134. 

[7] R. Liptser and V. Spokoiny, On Estimating a Dynamic Function of a 
Stochastic System with Averaging, Statistical Inference for Stochastic 
Processes 3(2000): 225-249. 

[8] W. Q. Zhu, Nonlinear stochastic dynamics and control in Hamiltonian 
formulation, ASME Appl. Mech. Rev. 59(4) (2006)230-248. 

[9] Y. Zeng, W. Q. Zhu, Stochastic Averageing of Quasi-Nonintegrable- 
Hamiltonian Systems Under Poisson White Noise Excitation, J. Apple. 
Mech.-Trans. ASME, 78(2011)021002-021011. 

[10] W. T. Jia, W. Q. Zhu, Stochastic averaging of quasi- 
non-integrable Hamiltonian systems under combined Gaus- 
sian and Poisson white noise excitations, Int J Nonlin Mech, 
http://dx.doi.Org/10.1016/j.ijnonlinmec.2012.12.003 

[II] Y. Xu, J. Duan and W. Xu, An averaging principle for stochastic dy- 
namical systems with Levy noise, Physica D. 240(2011), 1395-1401. 

[12] A. N. Kolmogorov, Wiener sche Spiralen und einige andere interessante 
Kurven im Hilbertschen, Raum, C.R.(Dokaldy) Acad.Sci.URSS(N.S.). 
26(1940), 115-118. 



22 



[13] B. B. Mandelbrot and J. W. Van Ness, Fractional Brownian motions, 
fractional noises and applications, SIAM Review 10(4)(1968), 422-427. 

[14] N. Chakravarti and K. L. Sebastian, Fractional Brownian motion mod- 
els for ploymers, Chemical Physics Letter. Vol. 267(1997), 9-13. 

[15] Y. Hu and B. 0ksendal, Fractional white noise calculus and application 
to finance, Infin.Dimens. Anal. Quantum Probab.Relat. Topics 6(2003), 1- 
32. 

[16] W. E. Leland, M. S. Taqqu, W. Willinger and D. V. Wilson,On the self- 
similar nature of ethernet traffic, IEEE/ACM Trans. Networking. 2(1994), 
1-15. 

[17] R. Scheffer and F. R. Maciel, The fractional Brownian motion as a 
model for an industrial airlift reactor, Chemical Engineering Science. 
vol.56(2001), 707-711. 

[18] I. Norros, E. Valkeila and J. Virtamo, An elementary approach to a Gir- 
sanov formula and other analytivcal resuls on fractional Brownian motion, 
Bernoulli. 5(1999),571-587. 

[19] A. N. Shiryaev, Essentials of Stochastic Finance : Facts , Models and 
Theory. World Scientific, New Jersey, 1999. 

[20] R. T. Baillie, Long memory processes and fractional integration in 
econometrics, Journal of Econometrics, 73(1996): 5-59. 

[21] G. Jumarie, Stochastic differential equations with fractional Brownian 
motion input. Int. J. Syst. Sci. 24(6)(1993), 1113-1132. 

[22] G. Jumarie, On the solution of the stochastic differential equation of 
exponential growth driven by fractional Brownian motion, Appl. Math. 
Lett. 18 (2005), 817-826. 

[23] S. C. Kou and X. S. Xie, Generalized Langevin Equation with Fractional 
Gaussian Noise : Subdiffusion within a Single Protein Molecule. Phys. 
Rev. Lett. 93(2004), 180603. 

[24] O. Y. Sliusarenko, V. Y. Gonchar, A. V. Chechkin,, I. M. Sokolov and R. 
Metaler, Kramers-like escape driven by fractional Brownian noise, Phys. 
Rev. E. 81(2010), 041119. 



23 



[25] W. Dai and C. C. Heyde, Ito formula with respect to fractional Brow- 
nian motion and its application, Journal of Appl. Math. and Stoch.Anal. 
9(1996), 439-448. 

[26] L. Decreusefond and A. S. Ustunel, Fractional Brownian motion : The- 
ory and applications, ESAIM:Proceedings 5(1998), 75-86. 

[27] D. Feyel and A. de la Pradelle, Fractional integrals and Brownian pro- 
cesses, Potential Analysis. 10(1996), 273-288. 

[28] P. Carmona and L. Coutin, Stochastic integration with respect to frac- 
tional Brownian motion, Ann. Inst. 39(2003), 27-68. 

[29] E. Alos and D. Nualart, Stochastic calculus with respect to the frac- 
tional Brownian motion, Ann. Probab. 29(2001), 766-801. 

[30] L. C. Young, An inequality of the Holder type connected with Stieltjes 
integratin, Acta Math. 67(1936), 251-282. 

[31] M. Zahle, Integration with respect to fractal functions and stochastic 
calculus II, Math. Nachr. 225(2001), 145-183. 

[32] F. Russo and P. Vallois, Forward, backward and symmetric stochastic 
integration, Probab. Theory Rel. Fields, 97:403-421. 

[33] F. Biagini, Y. Hu, B. Oksendal and T. Zhang, Stochastic Calculus for 
Fractional Brownian Motion and Applications, Springer- Verlag, London, 
2008. 

[34] T. E. Duncan, Y. Hu and B. Pasik-Duncan, Stochastic calculus for 
fractional Brownian motion I: Theory. SIAM J. Control Optim. 38(2000) 
582-612. 

[35] Y. S. Mishura, Stochastic Calculus for Fractional Brownian Motion and 
Related Processes, Springer- Verlag, Berlin, 2008. 

[36] T. Lyons, Differential equations driven by rough signals, Rev. Mat. 
Iberoamericana. 14(1998), 215-310. 

[37] D. Nualart and A. Rascanu, Differential equations driven by fractional 
Brownian motion, Collect. Math. 53(2002), 55-81. 



24 



